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Abstract. A nonlinear parabolic differential equation with a quadratic non- 
linearity is presented which has at least one equilibrium. The linearization 
about this equilibrium is asymptotically stable, but by using a technique in- 
spired by H. Fujita, we show that the equilibrium is unstable in the nonlinear 
setting. The perturbations used have the property that they are small in every 
IP norm, yet they result in solutions which fail to be global. 



1. Introduction 

This article demonstrates that in infinite-dimensional settings, stability of the 
linearization about an equilibrium of a dynamical system is not sufficient to ensure 
that the equilibrium is stable. This is in stark contrast to the situation in finite- 
dimensional settings, where stability of the linearized system implies stability of 
the equilibrium. (See [T], for instance.) A crucial point is that the system exhibited 
has a spectrum which includes zero, so stability is possible (as in the unforced heat 
equation), though not guaranteed. 

We study classical solutions to the Cauchy problem 

f JL«(t, a,) = £^ u ( t , x) - 2f{x)u(t, x) - u 2 (t, x) 

(1) I u(0,x) = h{x) € C°°(R) 
|i>0,iel, 

where / € Cq°(R) is a positive function with two bounded derivatives. (By Cq°, we 
mean the space of smooth functions which decay to zero.) Since the linear portion 
of the right side of dJ) is a sectorial operator, we can use it to define a nonlinear 
semigroup. [6] [10] The standard regularity theory for parabolic equations turns 
([l} into a smooth dynamical system, the behavior of which is largely controlled 
by its equilibria. This problem evidently has at least one equilibrium, namely 
u(t, x) = for all t, x. Depending on the exact choice of /, there may be other 
equilibria, however they will not concern us here. The linearized form of (TT]) about 
this equilibrium is evidently 

(2) ^u{t,x) = -^u{t,x)-2f{x)u{t,x). 

The zero function is asymptotically stable for the linearized problem, by a standard 
comparison principle argument. [3j However, using a technique pioneered by Fujita 
in [5], we will show that this equilibrium is not stable in the nonlinear problem, 
even if the initial condition has small p-norm for every 1 < p < oo. Fujita showed 
that if / = 0, then the zero function is an unstable equilibrium of {!]). The cause 
of the instability in ([T]) is the decay of /, for if / = const > 0, then the comparison 
principle shows that the zero function is stable. We extend Fujita's result, so that 
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roughly speaking, since / — ► away from the origin, the system is less stable to 
perturbations away from the origin. Another indication that there may be insta- 
bility lurking (though not conclusive proof) is that the decay of / means that the 
spectrum of the linearized operator on the right side of includes zero. [9] 

2. Motivation 

The problem |T]) describes a reaction-diffusion equation [3], or a diffusive logis- 
tic population model with a spatially-varying carrying capacity. The choice of / 
positive means that the equilibrium u = describes a population saturated at its 
carrying capacity. Without the diffusion term, this situation is well known to be 
stable. The decay condition on / means that the carrying capacity diminishes away 
from the origin. 

The spatial inhomogeneity of / makes the analysis of ([lj much more compli- 
cated than that of typical reaction-diffusion equations. The existence of additional 
equilibria for (JTJ is a fairly difficult problem, which depends delicately on /. (See 
[2] for a proof of existence of equilibria in a related setting.) 

3. Instability of the equilibrium 

Given an e > 0, we will construct an initial condition h € C°°(M) for the problem 
(JTJ), with \\h\\ p < e for each 1 < p < oo, such that ||u(f)||oo — > oo as t — -> T < oo. 
In particular, this implies that u = is not a stable equilibrium of ([1]), at least 
insofar as classical solutions are concerned. We employ a technique of Fujita, which 
provides sufficient conditions for equations like |T]) to blow up. [5] (Additionally, [3] 
contains a more elementary discussion of the technique with a similar construction.) 
Our choice for h can be thought of as a sequence of progressively shifted gaussians, 
and we will demonstrate that though each has smaller p-norm than the previous, 
the solution started at h still blows up. 

3.1. The technique of Fujita. The technique of Fujita examines the blow-up 
behavior of nonlinear parabolic equations by treating them as ordinary differential 
equations on a Hilbert space. Suppose u(t) solves 

du(t) 

(3) ^m =Lu (t )+N ( u ( t it), 

where L is a linear operator not involving t, and N may be nonlinear and may 
depend on t. Suppose that v(t) solves 

(4) 

where L* is the adjoint of L. Let J(t) = (v(t),u(i)). We observe that if \J(t)\ — > oo 
then either ||u(t)|| or ||u(i)|| also does. So if v(t) does not blow up, then we can 
show that ||u(t)|| blows up, and perhaps more is true. If we differentiate J(t), we 
obtain the identity 

pit) - f t (v(t),u(t)) 

I dv , A / , . du 

= (-L*v(t),u(t)) + (v(t),Lu(t) + N(u(t),t)) 
= («(i),A( U (<),i)), 
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where there is typically a technical justification required for the second equality. 
It is often possible to find a bound for (v(t), N(u(t),t)) in terms of J(t). So then 
the method provides a fence (in the sense of [7]) for J(t), which we can solve to 
give a bound on |J(t)|. As a result, the blow-up behavior of u{t) is controlled by 
the solution of an ordinary differential equation (for J(t)) and a linear parabolic 
equation (for v(t)), both of which are much easier to examine than the original 
nonlinear parabolic equation. 

3.2. Instability in LP for 1 < p < oo. We begin our application of the method 
of Fujita by working with L = 4^ — 2f and N(u) = —u 2 in ©. Since Q is 
then not well-posed for all t, we must be a little more careful than the method 
initially suggests. For this reason, we consider a family of solutions v c to (|4]) that 
have slightly extended domains of definition. It will also be important, for technical 
reasons, to enforce the assumption that the first and second derivatives of / are 
bounded. 

Definition 1. Suppose w — w(t,x) solves 



Define v e (s, x) = w(t — s + e, x) for fixed t > and s < t + e. Notice that by the 
comparison principle, v e (s,x) > 0. 

Lemma 2. Suppose that w solves ([5]). Then w, ^ S Co(R). 

Proof. The standard existence and regularity theorems for linear parabolic equa- 
tions (see [10., for example) give that w, §|, ^ S L 2 (K) and that w G C* 2 (JR). The 

comparison principle, applied to and Jj§^t gives that the first and second 

derivatives of w are bounded for each fixed t. (This uses our assumption that / has 
two bounded derivatives.) 

The lemma follows from a more general result: if g S C 1 n L P (M.) for 1 < p < 00 
and g' € L°°(R), then g £ Co(K). To show this, we suppose the contrary, that 
lima^oo g(x) (and possibly doesn't exist). By definition, this implies that 
there is an e > such that for all x > 0, there is a y satisfying y > x and \g(y)\ > e. 
Let S = {y\ \g{y)\ > e}, which is a union of open intervals, is of finite measure, and 
has supS* = 00. Let T = {y\ \g(y)\ > e/2}. Note that T contains S, but since g' is 
bounded, for each x € S, there is a neighborhood of x contained in T of measure 
at least e/Ug'Hoo- Hence, since supT = sup 5* = 00, T cannot be of finite measure, 
which contradicts the fact that g S L P (M.) with 1 < p < 00. □ 



Lemma 3. Suppose u : [0, T) x M — > R is a classical solution to {1} with u < and 
u(t) e L°°(R) for each t £ [0, T). Then 




(G) 




where w is defined as in Definition^ 



Proof. Define 
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First of all, we observe that since u € L°°(R), v e (s, -)u(s, •) is in for each 

s < t. 

Now suppose we have a sequence {m n } of compactly supported smooth functions 
with the following properties: [8] 

. m„ e C°°(R), 

• m n (x) > for all x, 

• supp(m„) is contained in the interval (— n — 1, n + 1), and 

• m n (x) = 1 for |x| < n. 

Then it follows that 



JAs) = lim / t> e (s, x)u(s, x)m n (x)dx. 

Now 

— J e (s) = — lim /" v e (s, x)u(s, x)m n {x)dx 
as as n— >oo J 

1 f 

= lim lim — / (vAs + h, x)u(s + h, x) — vJs, x)u(s, x)) m n (x)dx. 

h—>0 n^oo h J 

We'd like to exchange limits using uniform convergence. To do this we show that 
1 r 

(8) lim lim — / (vAs + h, x)u(s + h, x) — f e (s, x)u(s, x)) m n (x)dx 

n^oo h^O h J 

exists and the inner limit is uniform. We show both together by a little computation, 
using uniform convergence and LDCT: 



lim lim — / (v e (s + h,x)u(s + h, x) — v e {s, x)u(s, x)) m n (x)dx 

n^oc h^O h J 

f ( d d \ 

= lim / [ —v e (s,x)u(s,x) + vAs,x)—u(s,x) ) m n (x)dx 

n^oo J \OS OS J 

= lu^J \--^v e (s,x) +2f(x)v e (s,x)j u(s,x)m n (x) + 
w e (s, x) I — — ku(s, x) — u 2 (s, x) ~ 2f(x)u(x) ) m n (x)dx 



Minkowski's inequality has that 

v e um n dx 
since v e , m n > 0. This gives that 



1/2 , . N 1/2 

< / v t \u\m n dx < ( / v e m n dx I ( / v e u 2 m n dx 
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< 



v € (s, x)u 2 (s, x)m n (x)dx 

(J v e um n dx) 2 
J v e m n dx 



(/ v e udx)' 



J v e midx ' 

hence the inner limit of 1(5)) is uniform. On the other hand, 

\v e (s,x)u 2 (s,x)m n (x)\ < u e (s,a;)||it(s)||^ G L\R) 
so the double limit of (|8|) exists by dominated convergence. Thus we have the fence 
(9) dJ e ( s ) < (J £ ( S )) 2 _ 

ds ~ ||w e (s)||i' 
We solve the fence © to obtain (note J e < 0) 

1 < _dJ e (s) 1 



|» e (a)||i " (J e (s)) 

1 11 

-ds < 



2 



|« e («)||i " J e (t) J e (0) 

1 , 1 

-ds < -- 



|» e («)||l " J e (0) 

Taking the limit as e — > of both sides of the inequality yields 



1 , V 1 f f l 1 



w{t,x)h{x)dx < / t—77 rrpds = / || , r]] ds 

\Ja / \Jo \\w{s)\\i 

as desired. □ 

Remark 4. Since we are interested in proving the instability of the zero function 
in ([Tjl. consider u(0,x) = h(x) = — e for e > 0. Then ([6|) takes on the simple form 

So in particular, || w(t) ||oo blows up if there exists a T > such that e f Q T l^-jiT^ ds > 
1. 

The stability of the zero function in JTJ) depends on the stability of the zero 
function in (O - the linearized problem. If the zero function in the linearized 
problem is very strongly attractive, say ||w(t)||i ~ e - *, then 

f = (1 - e"*) < 1, 

Jo e 

and so a small choice of e < 1 does not cause blow-up via a violation of (fTU)) . On 
the other hand, blow-up occurs if it is less attractive, say ||u;(f)||i ~ t~ a for a > 0. 
Because then 

* s a t 

— ds = , 

t a a + 1 

whence blow-up occurs before t = 2±i, 
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In the particular case of f(x) = for all x, we note that w is simply a solution 
to the heat equation, which has = ||wo||i for all t (by direct computation 

using the fundamental solution, say), so blow up occurs. Thus we can recover a 
special case of the original blow-up result of Fujita in [5]. 

Theorem 5. Suppose a sufficiently small e > is given. Then for a certain choice 
of initial condition h(x) with \\h\\ p < e for all 1 < p < oo, there exists a T > for 
which lim t ^ T - ||u(t)||oo = oo. 

Proof. First, it suffices to choose ||u(0)||i < e and ||u(0)||oo < since 

i/p 

, 

< e. 



(| kW) "<IMir i)/p IMIi /p 



(11) — ^ = K, 



We assume, contrary to what is to be proven, that ||u(i)||oo does not blow up for 
any finite t. In other words, assume that u : [0, oo) x M — > R is a classical solution 
to (fTJ), with ||u(t)||oo < oo for all t. We make several definitions: 

• Choose < P < min |e, yj^- j- 

• Choose 7 > small enough so that 

A 

27 7 2 

for some some arbitrary K > 1. 

• Since < / € C§°(M.), we can choose an x\ such that 

(12) f( x ) < 7 when x < x\. 

• Next, we choose xq < x\ so that 

(13) Vi||/|u(l-erf(^E^ < 7 

for all < t < ^iy. Notice that any choice less than Xq will also work. 

• Choose the initial condition for (fTJ) to be 

(14) u(Q, x) = h(x) = -f}e p3 ' 2{x - x ^ 2 . 

This choice of initial condition has ||u(0)||oo = P < e, ||u(0)||i = 2n 1 / 2 f3 1 / i < 



e, and 



d 2 u(0) 



= /J, = 2/3 5 / 2 . (The value of /i will be important shortly.) 

• Finally, let wo(y) = S(y — xq) (the Dirac (5-distribution), and suppose that 
w solves ([5]). In other words, choose w to be the fundamental solution to 
([5| concentrated at x$. Note that the maximum principle ensures both that 
w(t,x) > for alH > and x G M and that \\w{t)\\i < ||w(0)||i = 1 for all 
t > 0. This allows us to rewrite ([6]) as 

(15) - t / w(t,x)h(x)dx < 1. 



Now we estimate the integral in (|15p . Notice that 



j t I w(t,x)(-h{x))dx = I -2f(x)w(t, x) )(-//(.<•) Ul.r 

~Q^2 + 2 f( x )H x ) ) w (^> X)dx, 
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where Lemma [2] eliminates the boundary terms. Now suppose z solves the heat 
equation with the same initial condition as w, namely 



(16) 



dz d z 

dt dx' 2 



z(0, x) — U>o(x) = S(x — Xq)- 



The comparison principle estabilishes that z(t, x) > w(t, x) for alH > and 
since f,w>0. As a result, we have that 



dt 



w(t,x) (—h(x)) dx > 



d 2 u 



dx 2 



+ 2f(x)h(x) z{t,x)dx 



> — [i — 2(3 I f(x)z(t,x)dx, 



where /i 
result, 



and (3 = ||u(0)||oo, which is an integrable equation. As a 



(17) / w(t, x) (-h(x)) dx > (3-fj,t-2(3 



o 

On the other hand using our choice for wq, 



1 (a;-y) 2 

f(x)—j^=e 4 = wo(y)dy dx ds. 
V Airs 



1 (.-»„)' 



f(x) p 1 4 S wo(y)dydxds— / / f(x) e dxds 

V Airs Jq J v Airs 



< 



[l I e 4=°' dx+ll/Hoo / e * *° } dx)ds 

V47TS 



4 2 

7 V< i 



< + A-erff^i^ 

2"! - Xq 



' ^ + ~*II/IU ( 1-erf. 

4 2 " J "°° V V 2Vi 



< 



we have used (fl"2"]) . (fH?|) . and assumed that < t < Then (fT5j) becomes 



1 > t / iu(i, x) (-h(x)) dx>/3t- [it 2 - 2(3jtVt = -2(3 5/2 t 2 - + (3t, 

\27K 



using our choices of fi, 7, and initial condition. Maple reports that the maximum 
of A(t) = -2(3 b/2 t 2 - 2 ^7K- /2 + is unic l ue ' occurs at < t a < , and has the 
asymptotic expansion 

A(t ) ~ K - !8Ky/(3 + A32K 3 (3 + 0{(3 Z/2 ). 

Thus for all small enough e > (3, we obtain a contradiction to (|15[) since K > 1. 
Thus, for some T < to < 00, lim t ^^- ||u(t)||oo = 00. □ 
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4. Discussion 

Theorem [5] gives a fairly strong instability result. No matter how small an initial 
condition to (TTJ) is chosen, even with all p-norms chosen small, solutions can blow 
up so quickly that they fail to exist for all t. This precludes any kind of stability for 
classical solutions. Like the analogous result in Fujita's paper, the kind of initial 
conditions which can be responsible for blow up are of the nicest kind imaginable 
- gaussians in either case! 

It must be understood that the argument in Theorem O depends crucially on 
the decay of /. Without it, the lower bound on J w(t,x)(—h(x))dx decreases too 
quickly. Indeed, if / = const > and h(x) > — /, then the comparison principle 
demonstrates that the zero function is asymptotically stable. On the other hand, 
any rate of decay for / satisfies the hypotheses of Theorem [SJ and so will cause ([T]) 
to exhibit instability. 

Finally, although we have examined the case where the nonlinearity in (fT]) is 
due to u 2 , there is no obstruction to extending the analysis to any nonlinearity 
like \u\ k , with degree k greater than 2. A higher-degree nonlinearity would result 
in a somewhat different form for ([6]), but this presents no further difficulties to 
the argument. Indeed, by analogy with Fujita's work, higher-degree nonlinearities 
would result in significantly faster blow-up. 
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